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We address the decomposition of a multi-mode pure Gaussian state with respect to a bi-partite 
division of the modes. For any such division the state can always be expressed as a product state 
involving entangled two-mode squeezed states and single mode local states at each side. The char- 
acter of entanglement of the state can therefore be understood modewise; that is, a given mode 
on one side is entangled with only one corresponding mode of the other, and therefore the total 
bi-partite entanglement is the sum of the modewise entanglement. This decomposition is generally 
not applicable to all mixed Gaussian states. However, the result can be extended to a special family 
of "isotropic" states, characterized by a phase space covariance matrix with a completely degenerate 
symplectic spectrum. 
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I. INTRODUCTION 

While the full characterization of bi-partite entangle- 
ment for mixed states is still an open problem, much is 
known for the case of pure states. Under a bi-partite 
division, any pure state may be written in the Schmidt 
form 



IAB 



^\<f>a)A\Xa)B 



(1) 



where (4> a \4>b) = (Xa\Xb) = $ab, with unique Schmidt 
coefficients y/p\. The bi-partite entanglement (entangle- 
ment entropy) can then be fixed uniquely by the asymp- 
totic yield [f| of maximally entangled states and becomes 
a function of the Schmidt coefficients only. Moreover, the 
Schmidt decomposition appears to have an "irreducible" 
structure: generally speaking, equation cannot be 
brought into a simpler form just by means of local trans- 
formations. For instance, a bi-partite system of n x n 
qubits cannot be generally brought to the form of a prod- 
uct of n entangled pairs under local unitary transforma- 
tions. 

However, in the context of Bosonic Channel Capacity, 
Holevo and Werner have shown that a multi-modeQ 
Gaussian mixed state can always be purified by enlarg- 
ing the system in such way that each normal mode is 
correlated with a corresponding single ancillary mode. 
This procedure achieves a pure Gaussian state between 
the system and ancilla in which the Schmidt decompo- 
sition takes the form of products of bi-partite two-mode 
Gaussian states. Implicit in these results is a general 
statement in the converse sense, which we believe is of 
considerable significance for the area of continuous vari- 
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able entanglement. The statement is that the bi-partite 
entanglement of multi-mode Gaussian pure states is in 
fact reducible to the product of entangled pairs of sin- 
gle modes. In other words, bi-partite entanglement of a 
Gaussian pure state is essentially I x 1 mode Gaussian 
entanglement. 



This result is directly applicable to various problems 
such as quantum-optical realizations of quantum infor- 
mation processing with Gaussian states 4] , and the char- 
acterization of the entanglement content of harmonic os- 
cillator chains[5j and bosonic quantum-fields 6, 7]. Con- 
sider for instance the vacuum state of a free scalar field, 
which is Gaussian. While one would expect that in this 
state, the structure of entanglement between a given re- 
gion of space and its complement would be of a rather 
complicated nature, such entanglement in fact occurs 
along separate "channels" , with each member of a set 
of collective modes in one region correlated with a corre- 
sponding single collective mode of the other. 



In this paper we present two different frameworks from 
which the modewise decomposition of Gaussian states 
can be deduced, and discuss some of its implications. In 
the following section, we present this modewise decom- 
position in the form of a theorem applicable to arbitrary 
entangled pure Gaussian states and show how it follows 
from properties of the Schmidt decomposition. In sec- 
tion lllll we deal with the case of mixed states. Using the 
correspondence between correlation matrices and Gaus- 
sian states, the modewise decomposition implies a cor- 
responding decomposition of covariance matrices. We 
therefore show how such a decomposition also holds for 
a certain class of "isotropic" Gaussian mixed states, de- 
fined from a corresponding symmetry of their covariance 
matrix. 
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II. MODEWISE DECOMPOSITION OF PURE 
GAUSSIAN STATES 

To begin with, suppose a collection of N canonical sys- 
tems or "modes" is partitioned into two sets, i.e., Alice's 
A = {A i} . . . , A m } and Bob's B — {Bi . . . , B n }, of sizes 
m and n respectively. If the quantum state of the modes 
is a pure Gaussian state \4>}ab, the following theorem 
characterizes the entanglement between Alice and Bob: 

Theorem 1: A Gaussian pure state \iP)ab f 0T rn + n 
modes A and B may always be written as 



I AB 



|Vl 



l^h B |0h JO). (2) 



for some s < min(m,n), where A = {A\ . . . , A m } and 
B = {B\, . . . ,B n } are new sets of modes obtained from 
A and B respectively through local linear canonical trans- 
formations, the states \^k) ire two-mode squeezed states 
fa/ of the form 



\1>i 



1 



(3) 



entangling the modes A/- and Bk for k < s, and |0)~ 

Ap 

and |0)^ are products of oscillator ground states for the 

remaining modes in A and B respectively. 

Before proving theorem 1, we first review some facts 
concerning the correspondence between Gaussian states 
and covariance matrices. Let us represent the canonical 
variables of a fc-mode system by the vector 



(4) 



where r\i is the two component vector r\i — (qi,Pi) T , and 
assume throughout that (rj) = for all states considered. 
A generally mixed Gaussian state p describing a system 
of fc-modes with (r}) — is completely specified by its 
covariance matrix (CM), defined as 



M = Re(n -q 1 



(5) 



A unitary transformation on p preserving the Gaussian 
character of the state implements a linear transforma- 
tion of the modes r\ = Sn, known as a symplectic (or 
linear canonical) transformation S S Sp(2k,W). Such a 
transformation preserves the canonical structure of the 
commutation relations [r],r] T ] = iJ 2 k, where the fc-mode 
symplectic matrix is given by 



■hk - J 2 , J 2 - ( _°i o 



and satisfies = — t-2k- Hence, a symplectic transfor- 
mation preserves the symplectic matrix under a similarity 
transformation, i.e., 

SJ 2 kS T = .hk => -(J2kS)(J 2k S T ) = t 2k - (6) 



Under such a transformation, the state p is brought to a 
new state p with CM M — SMS T . In particular, there 
exist symplectic transformations bringing the CM to the 
so-called Williamson normal form (WNF) 0, 



W = Ail 2 © A2II2 



Xkt 



k.U-2 



(7) 



where Xi are the non-negative eigenvalues of the matrix 
iJkM, also known as the symplectic eigenvalues. Ex- 
pressed in the product Hilbert space corresponding to 
the new set of modes {rji} for which W = Re(ry^ T ), the 
Gaussian state p acquires the particularly simple form 



p = pi ® p 2 ® ... ® pk 



(8) 



where pi is an oscillator thermal state for the i-th mode 

e-ft^ _ 1 
Tr(e-A^) ~ 



Pi 



^e-^ n \n)i{n\i. (9) 



Here, N± = afa,i is the number operator associated with 
ai = (qi + ipi)/y/2, and /3j is related to the symplectic 
eigenvalue Xi by & = ln[(A 4 + l/2)/(A 4 - 1/2)]. Note that 
as a consequence of the uncertainty principle, admissible 
Gaussian states satisfy the condition Vi , Aj > \ , with 
pure Gaussian states when Vz,Aj = \. For A,; = 1/2, 
Pi = |0)j(0|j is obtained as the limit of © as /3; — > 00. 

We now proceed with the proof of theorem 1. The 
Schmidt decomposition (J) automatically yields the di- 
agonal form of the partial density matrices for A and B: 



PA=^Pa\<Pa){<j)a\ PB = ^ Pa \Xa) (Xa \ 



(10) 



which are seen to be of equal rank and spectrum thus 
showing that the p a s are unique. The basis states \(j> a }A 
and \xo)b are also unique (up to phase factors) for non- 
degenerate p a and otherwise may be chosen to be ele- 
ments of any orthonormal basis spanning the degenerate 
subspace. Now, if \iP)ab is Gaussian, then the reduced 
density matrices are also Gaussian. Thus, pa and pb 
can be written in the form (JSJ) in terms of the set of 
modes bringing the local covariance matrices into WNF. 
Suppose that there are s modes in A and t modes in B 
with symplectic eigenvalue A ^ 1/2. Since the remaining 
modes factor out from the respective density matrices as 
projection operators onto their ground state, we may fac- 
tor \iP)ab as $) ab\ 0), |0) where 1 0) and |0) are 
collective ground states onto the modes with A = i and 
\ip) ab is the generally entangled state for the remaining 
modes, Ai...A s and B\, ..B t . Concentrate then on \i/j)ab, 
the partial density matrices of which may be written as 



PA = ? , ^ ry(A) \nA){n A \ 



51 Z (A) 
n A 



PB 



-pB-nB 

B) i n B )(n B \ 



Zi 1 



where ua — {n A ,...,n A } T and Hb = {n^ ,...,n s } T 
are s and t-dimcnsional vectors representing occupa- 
tion number distributions on each side and (3a = 



3 



{P Al ,~'Pji a } T and @ B = {\'-'^ t } T represent the 
distribution of thermal parameters on each side. Now, 
by our previous discussion, both density matrices have 
the same rank and the same eigenvalues. This means 
that there must exist a one-to-one pairing between the 
occupation number distributions Ha and Hb, and such 
that 



Pa ■ rt A = (3 B ■ n B . 



(11) 



We now observe that the pairing ft a is a homoge- 

neous linear map, since Ha — and Hb = are paired 
(all P's ^ 0) and {Ha + H' A ,H B + H' B ) satisfies <jTT|) if 
(Ha,Hb) and (H' A ,H' B ) satisfy However, if a lin- 

ear map is one-to-one then the domain and range have 
the same dimensions. Thus we see that s = t, in other 
words, the number of modes in A and B with symplectic 
eigenvalue different from 1/2 are the same. Now, label 
the modes on each side in ascending order of P, so that 
0<(3 Ai <0 A2 < ... </3 S2 < ...<0 Bs . 

Consider first the case Ha — {1,0, ..0} T , yielding the 
smallest non-zero value of (3 ■ Ha- By construction, this 
distribution must be paired with the smallest non-zero 
value of Pb • Hb, which is (or can be taken to be in the 
case of degenerate /3~ ) Hb — {1,0, ..0} T . We thus find 
that (|11|) has a solution provided that P A = (3~ (hence 
A _ = A- ), and by the linearity property we find for any 
Ha, the map n A — » n~ = n A . At this point we can 
repeat the procedure but applied to the subspace of the 
remaining modes, in other words, solve for a map be- 
tween H' A = {0,n A2 ,..,n A J and H' B = {0, , .., } 

such that /3a ■ H' A = Pb ■ n' B . By a similar argument we 
find that 13- = (3~ and n~ = n_ . Iterating the pro- 

' A 2 r B 2 A 2 B 2 b t> 

cedure until all the components are exhausted, we find 
that the admissible solutions to 111|) are Ha — Hb (with 
a freedom of re-ordering the labels of degenerate modes), 
provided that Pa = Pb- Reconstructing the Schmidt 
decomposition of |V>)ab from pa and ps we see that 



I ab — 



1 v - _l 



E 



|n>- \n). 

Z i i A . i i i 



(12) 



Thus, \iP)ab = Wab\0) a \0) n is of the form (2). 



III. ISOTROPIC GAUSSIAN MIXED STATES 

Although the previous result may be proved directly 
from general features of the Schmidt decomposition, it 
may also be embedded in the more general framework of 
mixed entangled Gaussian states. Theorem 1 can thus 
be seen to be a special case of a more general modewise 
decomposition theorem for a certain family of Gaussian 



mixed states. To define this family, we shall say that 
a covariance matrix is isotropic if there exists a sym- 
plectic transformation of the modes W — SMS T with 
S £ Sp(2k, M) that brings M to the form 



W = X Q 1 



2k 



An > - 



(13) 



An isotropic Gaussian state may thus be defined as a 
Gaussian state with an isotropic CM (An example of such 
a state would be the thermal state of a set of oscillator 
modes with degenerate frequencies). Note that all pure 
Gaussian states are isotropic with An = \ (H = 1). 

The more general theorem is a consequence of the fol- 
lowing Lemma concerning isotropic CMs: 

Lemma 1: Let M be an isotropic CM for m + n modes 
i] = rjA © r\B with symplectic eigenvalue An. Then there 
exist local symplectic transformations Ha — Sa^a o,nd 
Hb = Sb^b such that upon appropriate pairing of the 
modes, the covariance matrix takes the form 



M = M 



M 



8 ... © M_ 



A 1 



2(n+m— s) 



(14) 

for some s < min(m, n), where M is an isotropic 
correlation matrix for the two-mode sector n~ ~ = r> - © 

J A i B i A i 

n - of the form 



M 




= A?-Ag. (15) 



The diagonal elements Xi in \15}) are at the same 
time the symplectic eigenvalues of the local CMs Ma = 
Re (t]a t/a) and Mb — Re (rjs r/g) differing from An, and 
the last block in \14\l gives the CM for the remaining 
modes. 

Given the correspondence between CMs and Gaussian 
mixed states, the following extension of Theorem 1 im- 
mediately follows: 

Theorem 2: An isotropic Gaussian state p AB of sym- 
plectic eigenvalue An for the m + n modes A and B may 
always be written in the form 



Pi 



(16) 



where the new modes {A{\ and {Bi} are obtained by lo- 
cal symplectic transformations from A and B, p A are 
mixed Gaussian two-mode states with CM of the form 
$Wj) , and and jo! ' are mixed states for the remain- 

Ap Bp 

ing modes in A and B respectively with diagonal isotropic 
CM of symplectic eigenvalue An . 

To prove Lemma 1, first perform local symplectic 
transformations Ha®Hb = (Sa © Sb)va © Vb bringing 
the local CMs Ma = Re(r?A^), Mb — Re{?]sr] B ) into 
WNF. The total CM thus obtained may be written as 



M = Reirjr] 1 



W A K 
K T W B 



(17) 



4 



with W A = and W B = ©? =1 A fii l 2 . We 

next note a useful fact regarding isotropic CMs. If H13|) 
is satisfied, then M — XqSS t for some symplectic trans- 
formation S. However, as is easily verified, S' = SS T is 
a symmetric symplectic transformation. Consequently, 
© implies that an isotropic fc-mode CM satisfies: 

~{JM) 2 = \ 2 l 2k . (18) 

Replacing Ijl7|l m t° GH)> an d using the fact that 
[W, J] — 0, the following equations are obtained: 

W 2 A - (J m k){J n K T ) = X 2 Q t 2m (19a) 
W%- (J n K T )(J m K) = A 2 1 2 „ (19b) 
-W A K + J m KJ n W B = 0. (19c) 

Consider then a 2 x 2 sub-block jfiTy = (r/- ri T ) oi K. 

i Bj 

Since and are diagonal, from equation (|19c(l we 
find that 

A, K,j A . (20) 

It is not hard to verify that unless A^ = A g , this equa- 
tion has no solution for other than = 0. Thus we 
find that modes in A and modes B with different sym- 
plectic eigenvalues are uncorrelated 

Next, let rj~ and r?~ stand for the modes in A and B 
with the same local symplectic eigenvalue A, and group 
the modes according to their eigenvalues so that M takes 
the Jordan form M = (J) A M\ where each block M x is 
the CM for the degenerate eigenmodes rj~ ®rj~. Con- 
centrating on a given A, assume that g A and g B are the 
degeneracies of A in the symplectic spectra of W A and 
W B respectively, so that M x may be written as 

Mir -I)- « 2i) 

Now note that M\ is also an isotropic CM with symplec- 
tic eigenvalue A Q . Substituting ll2*TJ) into fH^I therefore 
yields 

■h gA KxJ 2gB Kj = (A 2 -A^)1 29A (22a) 
■h gB KlJ 2gA K x = (A 2 -A 2 )1 29B (22b) 
J2 9A K x J 2gB = K x . (22c) 

Taking the trace of equations (|22a|l and (|22b|l and using 
the cyclic property of the trace TT[J 2gA K\J 2gB Kj] = 
Ti[J 2gB KlJ 2gA K x }, we obtain that (A 2 - \ 2 ){g A - g B ) = 
. Thus we see that the respective degeneracies of the 
symplectic eigenvalue A in the local covariance matrices 
M A and M B must be the same for A 7^ Ao- In such a 
case, let g A = g B = g and from equations (|22|l deduce 
that 

KlJ 2g k A = K x J 2g kl = -(X 2 -X 2 )J 2g (23a) 
J2 a k x J 2g kl = K X K{. (23b) 



Next, define a matrix j3 = ©^03, where 03 is the standard 
Pauli matrix, and hence satisfying (3 2 = l 2g , (3 T = (3, 
J2g} = 0, and !3J 2g j3 — —J 2g - Re-expressing K in 
terms of some other 2g x 2g matrix 0\ as 

K x = ^X 2 -\ 2 O x (3 (24) 

and substituting into (|23|l we find that O x must satisfy 

O x J 2g O T x = J 2g , T x O x = l 2g ; (25) 

in other words, O x must be an orthogonal symplectic 
transformation. One can then perform a one-sided sym- 
plectic transformation, say on the A-modes, leaving 

the local CMs invariant in M x and bringing K x into the 
block diagonal form, i.e., 

K' x = T X K X = ^A 2 -A 2 a 3 . (26) 

i 

In this way, we achieve a pair-wise decomposition of the 
degenerate subspace itself, where each pair has a covari- 
ance matrix of the form (8). Finally, we note that for the 
degenerate subspace associated with A = A Q , in which 
the degeneracies on each side are not restricted to be the 
same, equations l|2^|) imply that k Xo k x = => k Xo — 
0. Therefore, local modes with symplectic eigenvalue A Q 
decouple, as expected from the pure case A Q = 5. 

To conclude with, we note that it is known that for 
a 1 x 1-mode Gaussian mixed state, the Peres-Horodecki 
partial transpose criterion|ll| is both necessary and suffi- 
cient ^3] for entanglement and hence for distillability^jj. 
In the present case of a two mode CM of the form Ijl5|l . 
the partial transpose criterion implies that the state is 
entangled iff Ai > Aq + \. Consequently, an isotropic 
Gaussian state is entangled and distillable iff at least one 
of the symplectic eigenvalues of its local CMs satisfies 
this condition. 
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